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Abstract 

n 

We construct the general solution of the equation wt+Yl w x k P^ k \w) = p(w) + [w, Tp(w)], 

k=l 

for the N x N matrix w, where T is any constant diagonal matrix, n, N £ N + and 
p( k \ p, p : M — > R are arbitrary analytic functions. Such a solution is based on the obser- 
vation that, as w evolves according to the above equation, the evolution of its spectrum 
decouples, and it is ruled by the scalar analogue of the above equation. Therefore the 
eigenvalues of w and suitably normalized eigenvectors are the iV 2 Riemann invariants. We 
also obtain, in the case p = p = 0, a system of N 2 non-differential equations character- 
izing such a general solution. We finally discuss reductions of the above matrix equation 
to systems of iV equations admitting, as Riemann invariants, the eigenvalues of w. The 
simplest example of such reductions is a particular case of the gas dynamics equations 

1 Introduction 



It is well-known that the general solution of the scalar first order quasi-linear Partial Differential 
Equation (PDE): 



+ J2 u ^P ik) ^ = P^' u:R n+1 ^R, (1) 



fc=i 



can be constructed by the method of characteristics [I] , converting (OQ) into the system of ODEs: 

S = p(«), (2) 

^ = p( k \u) ,k = l,..,n, [A) 

defined on the characteristic curves described by equations ([2b). If, in particular, p = 0, then 
the characteristic curves become straight lines and equations ([2]) are integrated in the form: 



x k = p( k \u)t + rj k , k = 1, .., n, 



where k = 1, .., N are arbitrary constants, / : M™ — > R is an arbitrary scalar function that 
can be identified with the initial condition: f{x\, .., x n ) = u(xi, .., x n , 0) of the Cauchy problem 
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in M. n , and u is constant along the characteristic straight lines. From equations ([3]) it also follows 
that the general solution u(x, t) of (pQ) for p = is defined implicitly by the non-differential 
equation 

u = f( Xl - p (1 \u)t, .., x n - p {n \u)t) . (4) 

Vector generalizations of equation p]) have been considered in several papers [21 [3j HI El El [7] . 
The main result of our work is the construction, using elementary spectral means, of the 
general solution of the following matrix generalization of equation (0Q): 

n 
k=l 

where w is the unknown NxN matrix function of the n+1 independent variables (cci, .., x n , t) G 
M™ +1 , T is any constant diagonal matrix, [•, •] is the usual commutator between matrices, and 
p( k \ p, p : R — ► R, k = 1, ..,n are n + 2 arbitrary scalar functions representable by positive 
power series, so that the quantities p^ h \w), k = 1, .., n, p(w) and p(u>) are well-defined functions 
of the matrix w. 

In §2 we show that, as w evolves according to (jHJ), the evolution of its spectrum decouples, 
and it is ruled by the scalar analogue of the above equation. Therefore the eigenvalues and 
suitably normalized eigenvectors of w are the N 2 Riemann invariants. In addition, the number 
of characteristic curves is only N (and not iV 2 , as it would be for a more generic vector system 
of iV 2 equations). 

In §3 we construct, in the case p = p = 0, the matrix generalization of the non-differential 
equation (j3J), given by a nonlinear system of iV 2 non-differential equations for the components 
of w. 

In §4 we introduce a distinguished reduction of equation ([5]) to a system of iV equations 
exhibiting, as Riemann invariants, the eigenvalues of w. 

We end this introduction with the obvious remark that solving equation (jSJ) is equivalent 
to solving the equation transposed to (jSJ) 

n 

™t + Yl P (k) (™)™*h = P^°) + [P(*) T > ™] > ( 6 ) 
fc=i 

whose solution is just the transposed of w: w = w T . 



2 The spectral solution of (El) 

The main characterization of the matrix equation ([5]) is that, although the iV 2 components 
w ij> hi = l)--)-^ °f w are coupled by the evolution, such a coupling disappears in spectral 
space, and the eigenvalues of w evolve according to the scalar analogue ([T]) of (J5J). 

Proposition 1 Let {e^\ .., e^} be the eigenvalues of matrix w, let {v^, .., v^} be the 
corresponding right eigenvectors, conveniently normalized to satisfy the iV conditions ifj = 

1, j = 1, .., N, where vf is the z-component of vector Then, if w evolves according to the 
matrix PDE the dynamics in spectral space is decoupled in the following way. 

i) Each eigenvalue evolves separately according to the scalar analogue of equation (J5j): 

n 

^' } + E e ^V J) )=P(e (j) ), J = l,-.,N. (7) 

k=l 
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ii) Known from equation (J7|), the corresponding normalized eigenvector vy> evolves accord- 
ing to the linear equation: 

n 

-E(^) P^^) = (Tj-T^Pi^), i,3 = 1,.,N. (8) 

k=l Xk 

Proof. The proof is by direct calculation. We substitute the spectral decomposition of w: 

w = vEv~ 1 , (9) 

into (jSJ), where E = diag(eW, .., e^^O) and v is the matrix of the normalized eigenvectors: 
Vij = v\ . After multiplying this equation from the left and from the right respectively by v^ 1 
and v, we obtain, after some manipulation: 



E t + Y,E x J k \E)-p{E) + 



fe=i 



v 



1 (v t + J2v ik) (E)+Tvp(E) ),E 



k=l 



(10) 



Separating the diagonal and off-diagonal parts of this equation, we obtain: 



E t +ZE x J k KE) = p(E), 

k n =1 (11) 
Vt+Yli v x k p^ k \E) + Tvp(E) = vh, h diagonal 
fe=i 

The diagonal matrix h is fixed by the normalization of matrix v. In our case, the diagonal part of 
v is the identity matrix; therefore, taking the diagonal part of flllb ). one infers that h = Tp(E). 
Equations (ITTh) and (TTTb) with h = Tp(E) are just equations d7j) and (|S]) respectively. □ 
The PDEs (J7J) and (jSJ) are converted into the system of ODEs: 

^ = p( e W), j = l,..,iV 

G) 



dv 



,„ -(Tj-TdvPfte®), i,j 1 V (12) 

** = p(*>( c W) ,fc = l,..,n, 

defined on the characteristic curves described by equations (|T2b). On such characteristic curves, 
the eigenvalues are obtained by the quadratures associated with (fT2k); once the eigenvalues are 
constructed, the eigenvectors are obtained solving the linear equations (|T2b). 
If p = p = 0, the characteristic curves become the straight lines 

x k = pW(e®)t + ri k , fc = l,..,7i, (13) 

where rjk are arbitrary integration constants, and each eigenvalue e^(xi, ..,x n ,t) is defined by 
the implicit equation 

e® = (rji, ..,Vn) = e 0) (x! - p (1) (e (j) )t, -,x n - p^(e^)t) , (14) 

where e^' : 1R™ — > R is an arbitrary function. For the normalized eigenvector constant 
along the characteristic curve (fTBT) . the general solution is explicit in terms of e^: 

yV\ Xl , ..,x n ,t) = v {j) (xi - p (1) (e (j) )t, - p {n \e^)t) , (15) 



where : W 1 — > JBL N is an arbitrary vector function. Then the general solution of the matrix 
system (JSJ) is achieved through the spectral formula ([9]). Such a solution depends, as it has to 
be, on the iV 2 arbitrary scalar functions j = 1,..,N and v\ , i,j = 1,..,N, i ^ j. We 
remark that, although the evolution in spectral space is decoupled, the components 
of w evolve along all the N characteristics, through the algebraic coupling given by 

If, in particular, one is interested in solving the Cauchy problem for the matrix system (J3J 
on M. n , for p = p = 0, one first remarks that, in this case, j = 1,..,N and v| , i,j = 
1, .., N, i 7^ j are identified with the initial conditions for the eigenvalues and eigenvectors: 



£ ^ {x\ , .., Xnj 6 (xi , . . , X n , 0) , 

v^xi, ..,x n ) = ..,x n ,0), 



(16) 



Therefore: i) given the initial condition 

f{xi,..,x n ) = w(x!,..,x n , 0), (17) 

one uniquely constructs the initial conditions e^>{x\, .., x n ) and v^^( ) for the eigenvalues 

and eigenvectors solving the associated eigenvalue problem: 

/v (i) = e w v y) , j = 1,..,JV (18) 

with the above normalization for the eigenvectors, ii) Given e^' and v^, one obtains the de- 
coupled evolution of the spectrum through the equations (fT4l) and (|T5|) . iii) The reconstruction 
of w(xi, ..,x n ,t) is finally achieved through the formula (Q. 

The elementary spectral solution of the nonlinear PDE ([3]) presented in this section is, to 
the best of our knowledge, new. 

The nonlinear PDEs ([5]) have been recently identified by the authors because they arise as 
the simplest examples of a class of nonlinear systems of PDEs in arbitrary dimensions generated 
by a novel dressing procedure based on a homogeneous integral equation with nontrivial kernel 



3 The analogue of the implicit equation (SI) 

Starting with the spectral solution of §5§ obtained in the previous section, it is possible to 
construct, if p — p — 0, the matrix analogue of the implicit non-differential equation (J3|). 

Proposition 2 Let fij : W 1 — ► IR, i,j = 1,..,N be iV 2 arbitrary scalar functions admitting 
formal positive power expansions, so that /^-(Mi, .., M n ) are well-defined NxN matrices, where 
Mi, .., M n are arbitrary matrices NxN. Denote by (fij(Mi, .., M n ))ki the (/^-component of the 
matrix /y(Mi, .., M n ). Then the general solution w of equation ([5]) is characterized implicitly 
by the following system of iV 2 non-differential equation: 

N 

Wij = - x n I - P [n \w)t)) k . ,1,3 = 1, .., N, (19) 

fc=i 

where Wij is the (ij')-component of matrix w and / is the NxN identity matrix. 
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Proof. Using equations ([9]), (J7J) and (IE]), we have that: 

= ElLiW^O* = Er =1 v*(x 1 -pW(eW)t,..,x n -p(")(eW)i)x 
e W( Xl _ p(i)( e «)t, . . . , Xn - p (n) (e (fc) )t)(v-%- { > 

Let / be the matrix having e^, I = 1, . . . , N as eigenvalues and v as matrix of eigenvectors, 



i.e.: fij = v^e^v^- 1 , then equation ( |20l) becomes: 



Wij 



1=1 

N 

v^(xx -p (1) (e (fc) )t, . . . ,x n -pW(e (/c) )t)^(^i -P [1 \e [k) % ■ ■ . ,x n - p^^t)^ 1 )^. 

k,l=l 

(21) 

Replacing now the definition of function of a matrix: 

/tf(x! - p«(e( fe ))t, ...,x n - P W(eW)t)5 ka = 

E (V" 1 )^ (/«(*!/ " P (1) (^K • ■ ■ ,Xn/ - ^WC V «' 
m,r=l 



into equation ( 1211) . we obtain the result. □ 

We end this section remarking that, as a byproduct of Propositions 1 and 2, the complicated 
system of N 2 non-differential equations (fT9l) defining the solution w is reduced to the solution 
of the single non-differential equation (fT4j) for the eigenvalues of w. 



4 A distinguished reduction 

In this section we briefly describe a reduction of the matrix equation (jSJ) exhibiting Riemann 
invariants coinciding with the eigenvalues j = 1, .., N. 

Consider the subspace of N x N matrices spanned by the basis {uq, ..,un-i} given by: 

= I, (wi)y = ■■ , {uk)ij = 8i+k,j, ■■ , (vN-i)ij = &i+N-ij, modiV. (23) 

This subspace is left invariant under matrix multiplication, since: 

UjUk = uJk^Jj = uJj+k, modiV, (24) 

therefore it defines a reduction of ([SD from the N 2 components of w to the N scalar coefficients 
Vk-, k = 1, .., N of the expansion: 

TV 

w = ^ VkUk-i- (25) 

k=l 

We remark that the mapping between the N dependent variables v^, k = 1, ..,N of the 
reduced system and the eigenvalues j = 1, .., N allows one to decouple completely the 
above dynamics; therefore the eigenvalues j = 1, .., N are the Riemann invariants for this 
reduced class of equations. 

In the remaining part of this section we show some examples of reduced systems. 

A system of 2 interacting fields in 2 + 1 dimensions is obtained, for instance, choosing 
N = 2, n = 2, pM(x) = x, p (2) (x) = ax 2 , p = p = 0: 

v\t + vivix! + v i v ix x + a{vf + v\)v\x 2 + 2av 1 v 2 V2 X2 = 0, , . 

V2t + v % v Xxi + ^1^2x1 + 2auxu 2 u lx2 + a{y\ + f|)f 2a;2 = 0. 
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If a = 0, equation ( 1261) reduces to a particular case of the well-known gas dynamics equations 
IB: 

v\ t + v x v Xxx + v 2 v 2xx = 0, ^ 

Vlt + V2Vl Xl + Vl"2 Xl =0. 

The eigenvalues of w: 

e« =v 1 + v 2 , e (2) =vi-V2, (28) 

evolve decoupled according to equation (j7j) and are the Riemann invariants of the above two 
systems (j2HP and (1271). 

The simplest example of a system of 3 interacting fields arises choosing N = 3, n = 
1, //^(x) = x, p = p = 0: 

u 2t + v 2 v lxx + i/!^ + u 3 u 3xi = 0, (29) 

^3t + V 3 U lxi + Wixx + "l^xi = 0- 

The three eigenvalues 

e (1) = vi + v 2 + z/ 3 , e (2) = | [2z/i - z/ 2 - ^3 + i\/3|f 2 - ^|] , ( ~n\ 
eW = \[2v l -v 2 -Vz-iV? > \p 2 -v z \] [M) 

evolve decoupled according to equation (J7J). But only is real, therefore the system is not 
hyperbolic and the eigenvalues and can be called Riemann invariants only in an extended 
sense. 

The systematic study of this reduced class and of other reductions of equation ([5]), in the 
search for integrable and applicative equations, is postponed to a subsequent paper. 
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